o 

(N 



Influence of Carrier-Carrier Scattering on Electron Transport in Monolayer 
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The influence of electron-electron scattering on the distribution function and transport characteristics of 
intrinsic monolayer graphene is investigated via an ensemble Monte Carlo simulation. Due to the linear 
dispersion relation in the vicinity of the Dirac points, it is found that pair-wise collisions in graphene do 
not conserve the ensemble average velocity in contrast to conventional semiconductors with parabolic energy 
bands. Numerical results indicate that electron-electron scattering can lead to a decrease in the low field 
mobility by more than 80 % for moderate electron densities. At high densities, the impact gradually diminishes 
due to increased degeneracy. 
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Since graphene was first realized experimentally in 
2004fi it has attracted significant interest due to its 
unique electronic properties. At low electron energies 
near the inequivalent Dirac points, there is no gap be- 
tween the valence and conduction bands and the dis- 
persion of the energy bands is linear.^ Extremely high 
electrical mobilities have been reported in suspended 
graphene, exceeding 10^ cm^/Vs near room tempera- 
turCj^ suggesting potential applications to ultrahigh fre- 
quency electronic devices.*^ The carrier density, which is 
controlled by the gate voltage, can be expected to vary 
by orders of magnitude. An interesting consequence of 
the linear energy dispersion is that the ensemble average 
velocity is not necessarily conserved upon an electron- 
electron (e-e) scattering event. Accordingly, inter-carrier 
collisions deserve careful consideration in the determina- 
tion of the transport properties. 

Das Sarma et al.^ found the inelastic e-e scattering rate 
and mean free path in graphene through the analysis of 
the quasiparticle self-energies. The scattering rate, cal- 
culated for electron densities from 1 — 10x10^^ cm~^ , was 
found to be of the same order of magnitude for electron- 
phonon scattering rates evaluated in the deformation po- 
tential approximation (DPA) i^ Several authors have con- 
sidered the electronic transport properties of graphene 
based on approaches such as the Monte Carlo simula- 
tion (see, for example, Ref. 0); however, the effects of e-e 
scattering have yet to be addressed to the best of our 
knowledge. In the present study, we examine the infiu- 
ence of this interaction mechanism in intrinsic monolayer 
graphene at room temperature. A full-band ensemble 
Monte Carlo method is used for accurate analysis of the 
distribution function as well as its macroscopic manifes- 
tations, particularly, the electron low-field mobilities and 
drift velocities. 

In both bulk and two-dimensional (2D) conventional 
semiconducting systems, e-e scattering has been well 
studied as is documented in the literaturcj^ During an 
e-e scattering event, both the energy and momentum 
are conserved. In a parabolic band structure, common 
in conventional semiconductors, momentum conservation 
directly leads to the conservation of velocity. This can 



be readily shown by multiplying the momentum conser- 
vation equation by h/ra, which gives vi + V2 = v'j^ + w'^^. 
It is then clear that, in a material with a parabolic band 
structure, e-e scattering has no direct effect on the drift 
velocity (which is an ensemble averaged quantity). The 
situation becomes different in graphene, where the energy 
dispersion in the region of the inequivalent Dirac points 
is linear, therefore removing the constraint that momen- 
tum conservation leads directly to velocity conservation. 
The problem of e-e scattering in nonparabolic bands 
was previously treated by Bonno et alM- Here we use a 
similar approach, making it specific to the linear band 
structure of graphene. The transition probability for two 
electrons, with wavevectors ki and k2, is given by Fermi's 
golden rule 
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s(ki,k2;k'i,k^)= — |Af|2^[i?(k'i) + i?(k^) 

-S(ki)-S(k2)], 



(1) 



where £'(k) is the electron energy dispersion relation and 
k'j^ and kj are the wavevectors of the final states af- 
ter scattering. The matrix element of interaction \M\, 
when accounting for exchange scattering (i.e., the indis- 
tinguishability of electron pairs with like spin after colli- 
sion), is given byi^ 



\M\ 



1 



[\V{q)\' + \Viq')\'-Viq)Viq')], (2) 



where the Coulomb scattering matrix V{q) between two 
electrons (ki — > k'^^; k2 — > 1^2) is 



V{q) 



27re^ 1 + cos(6ikik', ) 1 + cos{9i,^k'^ ) 
e{q)qA 2 2 



(3) 



with q = |ki — k'j^l, e{q) is the static dielectric function, 
^kk' is the angle between the wavevectors k and k', and 
A is the normalization area. The corresponding element 
V{q') (ki — ;> k2; k2 — > k'j^) can be expressed likewise with 
ka I . The dielectric function, which is valid for 



q' = |ki 

densities larger than n^lxlO^^ cm 
the random-phase approximation^^ 



is considered in 




FIG. 1. Schematic representation of the initial (ki,k2) 
and final (k'i,k2) electron pair states which satisfy energy 
and momentum conservation in the linear band structure of 
graphene. /? is defined as the angle between the long axis and 
the line OC and /?' between the long axis and the line OC 
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FIG. 2. (Color online) Electron-electron scattering rate 
^^, for electron concentrations n — IxfO^^ cm~^ (circle), n = 

cm~^ (square), and n = f xfO^"^ cm~^ (triangle). 



The scattering rate, including the effect of degeneracy, 
is then obtained as the double sum over all the final states 
k'j^ and the partner electrons k2 (with kg determined by 
momentum conservation) 
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Tee(ki 



=EE*(ki,k2;k'i,k^)/(k2) 

k2 k; 

x[i-/(k;)][i-/(k^)]. 



(4) 



where the /'s are the occupation probabilities. Within 
the Monte Carlo simulation, the Pauli exclusion principle 
is treated using the standard rejection technique. 

In graphene, the energy dispersion of the conduction 
band is linear in the vicinity of the Dirac points; i.e., 
E — hvf\k.\, with the Fermi velocity Vf (~ 10* cm/s). 
Therefore, the electrons behave as massless particles and 
can not be analyzed by the center of mass concept J^ 
However, by considering the equations of momentum and 
energy conservation, it can be shown that all the possible 
final states after scattering lie on an ellipse defined by 



ki + ka = k'l 



\i!-2 and 



ki| + |k2| = Ik'il 



As 



shown in Fig. [l] the initial point (i.e., the coordinate 
centrum at the K or K' point in the Brillouin zone) and 
the terminal point of the vector ki + ka specify the two 
foci of the ellipse, with the length of long axis of the 
ellipse determined by |ki| + jkaj. 

Our Monte Carlo model applied in the simulation in- 
cludes the full electronic band structure as calculated in 
the tight binding approximation. The electron-phonon 
scattering rates and phonon dispersion relations, previ- 
ously obtained from a first principles approach based on 
density functional perturbation theory, are included for 
all six branchesJ^ All electron-phonon transition pos- 
sibilities, including intravalley and intervalley, are ac- 
counted for as well as the e-e scattering. Finally, since our 
simulations include large electron densities we account 
for degeneracy, where Pauli exclusion is non- negligible, 
by implementation of the rejection technique in the se- 
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FIG. 3. (Color online) Drift velocity for n = IxfO^^ cm"^, 
n = 5x10^^ cm"^, and n = IxIO" cm"^ with (dashed line) 
and without (solid line) e-e scattering. 



lection of the final state after scattering»i^ 

In Fig. [21 we show the calculated e-e scattering rate, 
assuming a Fermi-Dirac distribution of initial states and 
that all final states are all available. The scattering rate is 
found to have a maximum of approximately 40-80 ps~^, 
for the energy range and electron densities considered, 
which is at least an order of magnitude larger than the 
total electron-phonon scattering rate.--- In Fig. [3l we ex- 
amine the effect of e-e scattering on the drift velocity 
from the low-field to saturation regimes. It is evident 
that there is a more than 80 % degradation of the low 
field mobility, which can be estimated by taking the slope 
of the drift velocity in the low field region, for an elec- 
tron concentration of n = 1 x 10^^ cm~^. Additionally, 
there is a degradation of the saturation drift velocity by 
^16 %. For low fields, when electrons rarely gain enough 
energy to emit an optical phonon, acoustic phonon and e- 
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FIG. 4. (Color online) Energy distribution function for n = 
1x10^^ cm~^ and n = 1x10^^ cm~^ with (dashed line) and 
without (solid line) e-e scattering. 



FIG. 5. (Color online) Calculated Edis/Etotai (solid line) and 
Vx/vf (dashed line) assuming the long axis of the ellipse is 
along the direction of the electric field. 



e scattering events dominate the electron dynamics and 
determine the mobility.^ As the field is raised, optical 
phonon scattering becomes more prevalent and the ef- 
fect of e-e scattering is somewhat mitigated. The dom- 
inance of e-e scattering at low energies, relative to the 
electron-phonon scattering rate, suggests that other scat- 
tering mechanisms (i.e., substrate dependent surface po- 
lar phonons, remote impurities, etc.) may tend to wash 
out the resulting degradation. 

For larger electron densities, n — 5 x 10^^ cm'^ and 
n = 1 X 10^^ cm~^, this effect is suppressed by the reduc- 
tion of available states with high and low energies. This 
can be seen from the electron energy distribution func- 
tions, given in Fig. 01 at an electric field oi E — 1 kV/cm 
for the highest and lowest carrier densities considered. 
Just as in conventional semiconductors, e-e scattering 
results in an extension of the thermal tail of the dis- 
tribution function, while simultaneously shifting its peak 
backward. Consequently, though an increase in electron 
density results in an increased e-e scattering rate, the dis- 
tribution function is less affected than in the moderate 
density case. This results in a larger degradation of the 
low-field mobility and saturation drift velocity for lower 
concentrations, as is evident in Fig. [3] 

Figure [5] presents how the selection of final state pairs 
effects the velocity of the ensemble which, as stated pre- 
viously, is not necessarily conserved. On average, the mo- 
mentum of an initial electron pair will be directed along 
the electric field which corresponds to the long axis of the 
ellipse, in Fig.[TJ The curve corresponding to i?diff/£^totai, 
where E^ie = {E'l—El^lis the energy difference of the final 
state electron pair, and i^totai — E[ + E'2, has maximal 
values when /3' is equal to or tt. Also shown in Fig. [SJ 
is the value of the final state drift velocity in the x di- 
rection, Vx — v[^ + v'2^, normalized by the fermi velocity 
Vf. Opposite to Editf / Etoted, the velocity in the direction 
of the electric field has minima at /3' is equal to or tt. 



This results in a decrease of the the ensemble drift veloc- 
ity for moderate electron densities, on the order of n = 
1x10^^ cm~^, since there are available high and low en- 
ergy states to scatter into. At larger densities, the lowest 
energy states are already mostly occupied, thereby pre- 
venting transition to the energy states which most greatly 
reduce the drift velocity. As a result, the effect of e-e 
scattering is most pronounced for n = 1x10"'^^ cm~^, and 
largely suppressed for n = 1x10^"^ cm^^. 

In summary, our results clearly indicate that due to e-e 
scattering, there is an electron density dependent degra- 
dation of the low-field mobility and saturation drift ve- 
locity in intrinsic monolayer graphene. This effect is most 
evident at moderate electron densities, while being sup- 
pressed at higher densities due to the unavailability of 
low energy states. For moderate electron densities, it is 
predicted that the low-field mobility and saturation drift 
velocity are degraded on the order of 80 % and 16 %, 
respectively. 
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